Observation of local chiral-symmetry breaking in globally centrosymmetric crystals. 
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A thorough tensor analysis of the Bragg-forbidden reflection (00.3)h in corundum systems having 
a global center of inversion, like V2O3 and a-Fe2 03, shows that anomalous x-ray resonant diffraction 
can access chiral properties related to the dipole-quadrupole (E1-E2) channel via an interference with 
the pure quadrupole-quadrupole (E2-E2) process. This is also confirmed by independent ab initio 
numerical simulations. In such a way it becomes possible, in this particular case, to estimate the 
intensity of the "twisted" trigonal crystal field (C3 symmetry) and, in general, to detect chiral 
quantities in systems where dichroic absorption techniques are ineffective. 



In recent years third generation radiation sources 
have made possible the detection of relatively small ef- 
fects in crystal electronic structure due either to mag- 
netic anisotropy^ or to interference between dipole and 
quadrupole (E1-E2) transition matrix elements. ^''^ For 
example, it is now well established that K-edge circu- 
lar dichroism in absorption is sensitive to the angular or- 
bital moment in magnetic spectroscopies like XMCD 
(x-ray magnetic circular dichroism), or to the peculiar 
physical quantity L^il^ in non-magnetic spectroscopies 
like XNCD-^'"* (x-ray natural circular dichroism). Here 
n = fAi — LAfis the toroidal (anapole^) orbital mo- 
ment, and z is the direction of the incoming photon. As 
apparent from the definition of il, the product L^^z is 
time-reversal even and inversion odd. Therefore if a para- 
magnetic system has a global inversion symmetry, even 
though at a local level this symmetry is broken, the re- 
sulting XNCD signal vanishes. A similar situation arises, 
mutatis mutandis, for XMCD in many antiferromagnets, 
where the locally broken time-reversal symmetry can be 
globally restored, making the total dichroic signal zero. 

A common way to circumvent this limitation is to use 
anomalous x-ray diffraction, where the local transition 
amplitudes are added with a phase factor that can com- 
pensate the vanishing effect due to the global symmetry. 
This technique has been widely used for antiferromagnets 
to study local magnetic effects.^ However to our knowl- 
edge no such technique has been reported to measure the 
quantity L^ilz for systems with a global inversion sym- 
metry. 

It is the purpose of this paper to show that this is pos- 
sible, by presenting a comprehensive theoretical analysis 
of the measured Bragg forbidden (00.3)?i "Finkelstein" 
reflection,^ in systems belonging to the corundum crystal 
class having a global inversion symmetry, namely V2 03^ 
and a-Fe203.^ In this latter compound, for example, this 
reflection has been interpreted^'^ as a pure E2-E2 pro- 
cess, giving rise to a six-fold periodicity in the azimuthal 
scan around the trigonal axis. We show here that such 
a description is not complete, as a non negligible E1-E2 
term is present, in V2O3 as well as in Q;-Fe203 and actu- 



ally in all the corundum systems. Its presence gives rise 
to a three-fold modulation of the azimuthal scan that is 
recognizable in the spectra of Paolasini et al.^, appears 
much more pronounced in Finkelstein et al.^ (see Figs. 2 
of both papers) and is related to the operator L^jJl^, as 
shown below. This modulation is also present in a more 
recent experiment by K. Ishida et al.^ on Q;-Fe203. 
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FIG. 1. Corundum and trigonal (rhombohedral) cells. 

We first focus our theoretical analysis on the paramag- 
netic corundum phase of V2O3 (space group R3c), so 
that we can neglect magnetic contributions. We shall 
comment on the general magnetic case later on. The 
four vanadium ions of the unit cell lie along the trig- 
onal axis (see Fig. 1), slightly offset from the centre 
of the distorted oxygen octahedra. Their positions in 
rhombohedral coordinates are^: Vi = {t,t,t), V2 = 
il-t,l-t,l- t), V3 = (1/2 + t, 1/2 + t, 1/2 + t), Vi ^ 
(1/2 -i, 1/2 -t, 1/2 -t), with t=0.1537. The four metallic 
sites are related by the following symmetry operations: 
V2=/Vi, V3=m3:Vi, V4=(72a;Vi, where / is the inver- 
sion in the reference frame of Fig. 1 and and rhx 
are, respectively, the two-fold rotation around a:-axis and 
the mirror plane orthogonal to it. 
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In the resonant regime the ii'-cdge elastic scattering 
amplitude, A(q,Lu), can be written in terms of the atomic 
scattering factors (ASF), fj{co), of the atoms at positions 
pj in the unit cell as: 



(1) 



where hq is the momentum transfer in the scattering 
process, Hlu the incoming and outgoing photon energy 
and the sum is over the four atoms in the rhombohedral 
primitive cell. The ASF is a second-order process in the 
electron-radiation interaction, whose explicit expression 
for core resonances reads, in atomic units: 



hoj - {En - Eq) - iVn 



(2) 



where I^'q"'"') is the ground state, with the origin taken on 
the j-th scattering atom, and Eq its energy; the sum is 
over all the excited states \^n), with corresponding en- 
ergies En- Finally r„ is a damping term that takes into 
account the core-hole and the finite life-time of the ex- 
cited states l^'ri) and the indices z, (,s) refer to the incident 
(scattered) properties of the photon field. Neglecting the 
electric dipole-magnetic dipole terms, which are usually 
small in the x-ray range, we c;an write the operator O up 
to the electric quadrupole contribution as: 



O = e • r(l — -k ■ r) 



(3) 



where e and k are the polarization and wave vector of the 
x-ray photon and r is the coordinate of the electron in 
the reference frame of the scattering atom. The matrix 
element in Eq. (2) depends only on the electronic part of 
the operator O, in such a way that the radiation param- 
eters, e and fc, can be factorised. After some algebra Eq. 
(2) can be written as a scalar product of two irreducible 
tensors'''^^: 



/,(^)=Sp,,(-)%(rt^b)(j-.^) 



(4) 



where Tq^^ depends on the incident and scattered polar- 
izations and wave vectors and Fq^^ is the tensor repre- 
senting the properties of the system. It is important to 



note that F!j^\j;uj) = {•^^q" IF^''' is such that 
it must belong to the totally symmetric representation 
of the local point group of the scattering atom (the Aig 
representation, in Bethe's notation). In fact it is easy to 

check that if i? is a symmetry operation for I^E'q"''') (ie. 



- /,t,W)if(p) 



{3)\ 
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In the corundum systems the local symmetry on each 
scattering atom is C3. Thus Eq. (5) imposes the fol- 
lowing restriction: only the irreducible tensors whose az- 
imuthal numbers with respect to the three-fold axis are 
or ±3 are allowed. A similar constraint was recognized 
for the case of haematite,^ although the analysis was 
there limited to an octahedral symmetry, loosing there- 
fore the E1-E2 contribution. If we choose the quantiza- 
tion axis coincident with the three-fold axis, the tensor 
components allowed by the C3 symmetry are: 



a) in the El-El channel: F^°\ F^^^ and i^i^^ 





(1) 







b) in the E1-E2 channel: F^'' , F^'^\ F^ , fP±F'_ 



p(3) 
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(3) 



c) in the E2-E2 channel: F^°\ F^i^K F^i^\ F^i^\ F^"*^ 



fP ± ^^3^ and FP ± F_ 



p(3) 



r(4) 



' ' ' 

(4) 

3 -L .t _3 aiiu i'3 -L i _3 

The tilded tensors refer to time-reversal odd (ie, mag- 
netic) quantities. As V2O3 in the corundum phase is 
paramagnetic, tilded tensors are all zero. Another ex- 
tinction rule comes from the structure factor. We are 
interested in the (00.3)/i reflection, (lll)r in the rhom- 
bohedral system. For (f=(lll)r, Eq. (1) becomes: 



/i + e-2-*/2 + e*"e2"*/3 



A{q,uj 

2cos(27rt)(l -m^)/i (El - El & E2 - E2) 
2i sin(27ri)(l - mx)/i (El - E2) 



-he''^e-2"*/4 = (e^"* + e-2''"/)(l - TOx)/i 



(6) 



where we uscicl the symmetry operations introduced 
above. The mirror rhx is such that (x,y,z)— >(-x,y,z). For 
the chosen reference frame its action on spherical tensors 
is given by-'^-^: rhxFq^^ = {—)^^^F^^, where P is the par- 
ity of the tensor (-|-1 for El-El and E2-E2 tensors and — 1 
for E1-E2 tensors). The combined action of glide-plane 
and Cs-symmetry forbids any El-El contributions, as 
expected, and leaves us with only three possible terms: 



A3) 



F^^, the one 



two come from the E1-E2 channel (Fg ' and F3 

and the other is of E2-E2 origin: F^'^^ 
recognized in Ref. 7. 

Their polarization and wavevector dependence shows 
that the signal is different from zero only in the (ttt 
channel and only for T^"^^ and Tg^"' - t!_^^. In fact the 

(3) (3) 

— T_3 term is proportional to k^ — k'^ or ky — k'y, 
where k and k' are the incident and scattered wave vec- 
tors and X and y are orthogonal to the trigonal axis. In 
the geometry of the (00. 3)?^, k^ = k'^ and ky = ky at any 
azimuthal angle. 

It is now easy to compare the theoretical azimuthal de- 
pendence around the momentum transfer Uq with the ex- 
perimental one^: the contribution of Tq^'* is constant with 
respect to the azimuthal angle (j) around the q-vector. 



{^q''\F^^\uj)\^q'') = {^q''\R~^F^^\ijj)R\^q^) (5) while that of T^'^^ — Tig'' is three-fold periodical in the 



.(4) 



Thus, the only allowed matrix elements are those of the 
components of F*-^-* (oj) that are invariant for any symme- 
try elements of the point group. 



amplitude.^ Note that time-reversal even quantities are 

real in the E2-E2 channel and imaginary in the E1-E2 
channel. Because of the imaginary unit in the E1-E2 
term of Eq. (6), both amplitudes are real and intefere. 
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Thus, the global dependence in the scattered intensity 
is proportional to (a + sin(3(^))^, ie, a three-fold modu- 
lation of the six-fold periodic intensity. The constant a 
takes into account the relative weight of the radial matrix 
elements for the two irreducible tensors. 

In order to estimate such a constant we have performed 
a numerical simulation based on the FDMNES package. 
The results are shown in Fig. 2 for the energy scan 
and in Fig. 3 for the azimuthal scan at the pre-edge 
energy (5465 eV). They both support our previous theo- 
retical considerations: the energy scan shows no features 
at El-El energies (as experimentally detected^) while a 
structure exists around the 3d-energies (5465 eV). Also 
the azimuthal scan of this structure, shown in Fig. 3 is 
compatible with the experimental data^ and makes clear 
its double-component origin. The possibility to evaluate 
separately the two contributions in the FDMNES pro- 
gram allows us to estimate the ratio between the constant 
E1-E2 signal and the maximum of E2-E2: a ~ 0.05. Note 
that this estimates includes the prefactors cos(27ri) and 
sin(27ri), for E2-E2 and E1-E2. For i = 0.1537 they are, 
respectively, 0.57 and 0.82. 
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FIG. 2. Numerical simulation for the energy scan of the 
(00.3)h reflection in corundum V2O3. Note the absence of 
any signal at the El-El 4p conduction band energies. Here 
and in the following figures the intensity is proportional to 
the square of the effective number of scattering electrons. 

A simple physical picture can help understanding the 
previous results. In the corundum system there are two 
main crystal fields that split the spherical degeneracy of 
the 3d electrons involved in the resonant process: the oc- 
tahedral crystal field is the biggest (lODq ~ 2 eV) and is 
responsible of the splitting to doubly degenerate e^ and 
triply degenerate X2g orbitals. The trigonal field, usu- 
ally one order of magnitude less, further splits the Xig or- 
bitals into doubly and singly degenerate levels. This field 
is made up of two contributions"'^'*: an axial field, of Q-^t, 
symmetry, describing the shift of vanadium ions along the 
trigonal axis, and a "twisted" field of C3 symmetry. This 
latter takes into account the small rotation of the oxy- 



gens around the trigonal axis, as well as the positions of 
the second-nearest-neighbors vanadium ions: both these 
latter features break the vertical mirror plane. In the ab- 
sence of a trigonal splitting only t'^'^ —T'^'I contributes to 
the signal, as it is totally symmetric in an octahedral po- 
tential. No inversion-breaking terms would be possible in 

this case. The detection of the term Tg oc L^^^ is made 
possible only by the reduction of the local symmetry to 
C3. Note that the allowed Tq ' tensor is time- reversal 
even and inversion-odd. The identification with L^^^ 
follows from the fact that this latter is the only phys- 
ical quantity detectable at the K-edge with these same 
properties. ^'"'^ Notice also that if the point symmetry of 
vanadium ions were C31,, as maintained by Carra et al.^'^ 
the term Lz^z would be forbidden by symmetry, as it 
changes sign under the action of a vertical mirror plane 
and no modulation could be seen in the azimuthal scan. 
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FIG. 3. Numerical simulation for the azimuthal scan 
around the (00.3)h direction at vanadium pre K-edge (5465 
eV). Also shown is the pure E2-E2 simulation in octahedral 
symmetry. In this last case the six-fold periodicity is clearly 
recovered. 

All this has been tested by a numerical simulation 
where a "straightened" corundum cell has been used, 
with a perfect octahedral environment around each vana- 

(2) 

dium ion. In such a condition, there is no Tq terms, as 
it should be. The E1-E2 signal, then, increases with the 
degree of "twisting" of the trigonal deformation of the 
oxygen octahedra. 

Similar results are obtained for the magnetic a-Fe2 03, 
as shown in Fig. 4. Neglecting for the moment mag- 
netic contributions, we observe that in this case the ratio 
between two consecutive maxima is rather pronounced 
(about 3/2), in keeping with the experimental findings.^ 
However it should be kept in mind that this ratio depends 
quite substantially on the photon energy, as numerically 
verified. If we now consider magnetic contributions, when 
the system is below the Morin temperature (Ta/ ~ 260 
K) the magnetic moments are parallel to the trigonal 
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axis and order in the sequency ti Ti ii J- respectively for 
ions 1,2,3,4 in Fig. 1. Thus the magnetic space group is 
again R3c and Eq. (6) still holds. Yet, time-reversal is 
locally broken, and this allows in principle contributions 
from the tilded tensors, both in the E1-E2 and E2-E2 
channels. From this point onward, tensor group analysis 
cannot be pushed further and we need quantitative calcu- 
lations. Our numerical simulations give a magnetic con- 
tribution about a factor 10~^ less than the non-magnetic 
one. This numerical result is in keeping with that of Ref. 
5, where the experimental data had been interpreted as 
non-magnetic, because the temperature dependence of 
the azimuthal periodicity was not affected by the tran- 
sition through the Morin temperature. In fact, if the 
three-fold signal were magnetic, it should change drasti- 
cally above T^v/, where all the spins align perpendicularly 
to the trigonal axis. In this latter case, while the three- 
fold lattice symmetry is still preserved, the magnetic one 
is broken, and more tilded tensors components, with dif- 
ferent periodicity, could in principle contribute. 
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FIG. 4. Same as Fig. 3 for Q-Fe203 (E=7105 eV). 

Unfortunately, the small ratio of the magnetic to non- 
magnatic signal makes difficult, at present, a systematic 
study of the magnetic space group based on this kind 
of analysis for other corundum systems like Cr203 or 
Ti203. Such a study would be interesting for many rea- 
sons: in the first case, in the light of the recent claim 
for magnetochirality,^^'^^ it could provide an indepen- 
dent verification. In the second case, because of an old 
controversy regarding the magnetoelectricity in Ti203,^* 
it could give more detailed information on the possible 
magnetic symmetries of the system. In all cases, the pres- 
ence/absence of tilded tensors would be determined by 
their different polarization and wavevector dependence, 
as well as by the changes in Eq. (6) due to the different 
magnetic space groups (in the case of Cr203, for exam- 
ple, both / and rhb appear with a minus sign for mag- 



netic tensors). In general for these and similar studies 
one should rely on better resolved spectra or on "pure" 
magnetic reflections. 

In conclusion, by a complete tensor analysis we have 
shown with a simple example that x-ray elastic resonant 
scattering and XNCD can be considered complementary 
tools for the detection of the chiral quantity Lzflz- In 
particular resonant diffraction, due to the possibility of 
adjusting phase factors, is the only technique that allows 
a detection of such a signal in systems where the inver- 
sion symmetry, though broken at a local level, is globally 
restored. 
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